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Abstract. This paper is concerned with some nonlinear propagation phe- 
nomena for reaction-advection-diffusion equations in a periodic framework. 
It deals with travelling wave solutions of the equation 

Mt = V ■ {A{z)\/u) + q{z) ■ Vu + f{z, u), t E R, z E n, 

propagating with a speed c. In the case of a "combustion" nonlinearity, the 
speed c exists and it is unique, while the front u is unique up to a translation 
in t. We give a min — max and a max — min formula for this speed c. On the 
other hand, in the case of a "ZFK" or a "KPP" nonlinearity, there exists 
a minimal speed of propagation c*. In this situation, we give a min — max 
formula for c*. Finally, we apply this min — max formula to prove a variational 
formula involving eigenvalue problems for the minimal speed c* in the "KPP" 
case. 
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1 Introduction and main results 

1.1 A description of the periodic framework 

The goal of this paper is to give some variational formula for the speeds 
of pulsating travelling fronts corresponding to reaction-diffusion-advection 
equations set in a heterogenous periodic framework. In fact, many works, 
such as Hamel [6], Heinze, Papanicolaou, Stevens [9], and Volpert, Volpert, 
Volpert [K)] treated this problem in simplified situations and under more 
strict assumptions. In this paper, we treat the problem in the most general 
periodic framework. We are concerned with equations of the type 

Ut = V ■ {A{z)Vu) + q{z) ■Vu+ f{z,u), t E R, z E n, 
V ■ A \/u{t, z) = 0, t E R, z E dn, 

where i>{z) is the unit outward normal on dQ at the point z. In this context, 
let us detail the mathematical description of the heterogeneous setting. 

Concerning the domain, let > 1 be the space dimension, and let d be an 
integer so that 1 < d < N. For an element z = (xi, X2, ■ ■ ■ , Xd, Xd+i, ■ ■ ■ , xn) G 
M^, we denote by x = (xi, X2, ■ ■ ■ , Xd) and hj y = {xd+i, ■ ■ ■ , xn) the two 
tuples so that z = (x, y). Let Li, ■ ■ ■ ,Ldhed positive real numbers, and let 
Q be a C'^ non empty connected open subset of satisfying 

{3R>0;y{x,y)En,\y\< R, 
V(A:i,--- ,kd) E L{L X ■■■ X L^Z, VL 

where {ei)i<i<N is the canonical basis of M^. 

Ks d > 1, one notes that the set Vt satisfying (1.2) is unbounded. We 
have many archetypes of such a domain. The case of the whole space 
corresponds to d = N, where Li, . . . , are any positive numbers. The case 
of the whole space with a periodic array of holes can also be considered. 
The case d = 1 corresponds to domains which have only one unbounded 
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dimension, namely infinite cylinders which may be straight or have oscillating 
periodic boundaries, and which may or may not have periodic perforations. 
The case 2 < d < N — 1 corresponds to infinite slabs. 

In this periodic situation, we give the following definitions: 

Definition 1.1 (Periodicity cell) The set 

C = {ix,y)e n; xi G (0, Li), ■ ■ ■ , G (0, L^)} 
is called the periodicity cell of Q. 

Definition 1.2 (L-periodic fields) A field w : ^ is said to be L- 
periodic with respect to x if 



w{xi + ki, 



+ kd,y) = w{xi, ■■■ ,Xd,y) 



almost everywhere in fl, and for all k = (fci, ■ ■ ■ , kd) G Y\ LiZ. 



i=l 



Let us now detail the assumptions concerning the coefficients in (1.1). 
First, the diffusion matrix A{x, y) = {Aij{x, y))i<i,j<N is a symmetric C^'^( fl ] 
(with S > 0) matrix field satisfying 



A is L— periodic with respect to x, 
30 < ai < a2;V(x,?/) G fi,V^ G M^, 



;i.3) 



^<i,j<N' 



The underlying advection q{x,y) = {qi{x,y), ■ ■ ■ ,qN{x,y)) is a C^'^{Q) 
(with 6 > 0) vector field satisfying: 



q is L— periodic with respect to x, 
V ■ g = in TI, 
q ■ u = on dQ , 



;i-4) 



VI < z < d, 



qi dx dy = 0. 
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Lastly, let / = /(x, y, u) be a function defined in x M such that 

/ is globally Lipschitz-continuous in 1] x M, 

V(x,2/)G n, VsG (-oo,0]U [l,+oo), /(s,x,2/) = 0, 
3 p G (0, 1), V(x, y) G n, V 1 - p < s < s' < 1, 
fix,y,s) > f{x,y,s'). 



(1.5) 



One assumes that 

/ is L— periodic with respect to x. (1-6) 

Moreover, the function / is assumed to be of one of the following two 
types: either 



3^ G (0,1), V(x,2/)G(^, VsG [0,^], /(x,2/,s) = 0, 
Vs G (6', 1), 3(x, y)Gr2 such that /(x, y, s) > 0, 
or 

3 5 > 0, the restriction of / to f2 x [0, 1] is of class C^''', 
Vs G (0, 1), 3 (x, y) G n such that /(x, y, s) > 0. 



;i-7) 
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Definitions 1.3 A nonlinearity f satisfying (1-5), (1-6) and (1-7) is called 
a "combustion" nonlinearity. The value 9 is called the ignition temperature. 
A nonlinearity f satisfying (1-5), (1.6), and (1.8) is called a "ZFK" (for 
Zeldovich- Frank- Kamenetskii) nonlinearity. 
If f is a "ZFK" nonlinearity that satisfies 

fl{x,y,0)= lim f{x,y,u)/u > 0, (1.9) 
«— >o+ 

with the additional assumption 

V(x,t/,s) G fix (0,1), < f{x,y,s) < fi{x,y,0) x s, (1.10) 

then f is called a "KPP" (for Kolmogorov, Petrovsky, and Piskunov, see [11]) 
nonlinearity. 

The simplest examples of "combustion" and "ZFK" nonlinearities are 
when f{x,y,u) = f{u) where: either 

/ is Lipschitz-continuous in M, 

3 ^ G (0, 1), /(s) = for all s G (-00, ^] U [1, +CX)), 
f{s) > for all s G (^,1), ^ ' 



3 p G (0, 1 — 6^), /is non-increasing on [1 — p, 1] 



or 



' / is defined on M, / = in M \ (0, 1), 
3 5 > 0, the restriction of / on the interval [0, 1] is C^''^([0, 1]), 
/(O) = /(I) = 0, and f{s) > for all s G (0, 1), ^ ' ' 

^ 3 p > 0, /is non-increasing on [1 — p, 1]. 
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If / = f{u) satisfies (1.11), then it is a homogeneous "combustion" non- 
linearity. On the other hand, a nonhnearity / = f{u) that satisfies (1.12) is 
homogeneous of the "ZFK" type. Moreover, a KPP homogeneous nonhnear- 
ity is a function / = f{u) that satisfies (1.12) with the additional assumption 

Vse(0,l), 0<f{s)<f'{0)s. (1.13) 

As typical examples of nonlinear heterogeneous sources satisfying (1.5- 
1.6) and either (1.7) or (1.8), one can consider the functions of the type 

f{x,y,u) = h{x,y)g{u), 

where his a globally Lipschitz-continuous, positive, bounded, and L— periodic 
with respect to x function defined in Q, and g is a. function satisfying either 
(1.11) or (1.12). 

Definition 1.4 (Pulsating fronts and speed of propagation) Let e = 

(e^, ■ ■ ■ ,e'^) he an arbitrarily given unit vector in M^. A function u = u{t, x, y) 
is called a pulsating travelling front propagating in the direction of —e with 
an effective speed c ^ 0, if u is a classical solution of: 

' Ut = V ■ {A{x, y)Vu) + q{x, y)-Vu+ f{x, y,u), t e R, (x, y) E n, 
V ■ A Vu{t, x,y) = 0, t e M, (x, y) G dQ, 

< yk G JjiviZ, y{t,x,y) G ]R X fi, u{t-\ ,x,y) = u{t,x + A;, y), (1.14) 

1=1 ^ 
lim u(t,x,y) = 0, and lim u{t,x,y) = 1, 

x-e— »— 00 x-e—>+oo 

< n < 1, 

where the above limits hold locally in t and uniformly in y and in the directions 
of M*^ which are orthogonal to e . 

Several works were concerned with pulsating travelling fronts in periodic 
media (see [1], [2], [10], [12], [14], [1.5], and [18]). 

In the general periodic framework, we recall two essential known results 
and then we move to our main results. 

Theorem 1.5 (Berestycki, Hamel [^]) Let VL be a domain satisfying (1.2), 
let e be any unit vector ofW^ and let f be a nonlinearity satisfying (1.5-1.6) 
and (1.7). Assume, furthermore, that A and q satisfy (1.3) and (I.4) respec- 
tively. Then, there exists a classical solution (c, m) of (1.14)- Moreover, the 
speed c is positive and unique while the function u = u{t,x,y) is increasing 
in t and it is unique up to a translation. Precisely, if {c^,u^) and (c^,m^) are 
two classical solutions of (I.I4), then = (? and there exists /i G M such 
that u^{t, X, y) = u^it + h, x, y) for all {t, x,y) E M x Q. 
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In a periodic framework, Theorem 1.5 yields the existence of a pulsating 
travelling front in the case of a "combustion" nonlinearity with an ignition 
temperature 6. It implies, also, the uniqueness of the speed and of the profile 
of u. For "ZFK" nonlinearities, we have 

Theorem 1.6 (Berestycki, Hamel [ ]) Let VL he a domain satisfying (1.2), 
let e he any unit vector in M'^ and let f he a nonlinearity satisfying (1.5-1.6) 
and (1.8). Assume, furthermore, that A and q satisfy (1.3) and (1.4) respec- 
tively. Then, there exists c^^^ j(e) > such that the prohlem (1.14) has no 
solution (c, m) such that ut > in M. x Q if c < c^,^,gj(e) while, for each 
c > CQ^qy(e), it has a solution {c,u) such that u is increasing in t. 

In fact, the existence and the monotonicity of a solution u* = u*{t, x, y) of 
(1.14) for c = CQ^A,q,fi^) > holds by approaching the "ZFK" nonlinearity / 
by a sequence of combustion nonlinearities (/e)^ such that fe^f uniformly 
in M X r2 as 6* \ O"*" (see more details in step 2 of the proof of formula 
(1.17) below, section 4). It follows, from Theorem 1.5, that for each ^ > 0, 
there exists a solution {c0,ug) of (1.14) with the nonlinearity fe such that 
ue is increasing with respect to t. From parabolic estimates, the functions 
ug, converge up to a subsequence, to a function u* in Cf^^{R x fl) as 9 
O"*". Moreover, Lemmas 6.1 and 6.2 in [1] yield the existence of a constant 
c*(e) = CQ^^ gj(e) > such that cg c*(e) as 6* \ 0. Hence, the couple 
(c*(e),M*) becomes a classical solution of (1.14) with the nonlinearity / and 
one gets that u* is nondecreasing with respect to t as a limit of the increasing 
functions uq. Finally, one applies the strong parabolic maximum principle 
and Hopf lemma to get that w is positive in M x In other words, u* is 
increasing in t G M. Actually, in the "ZFK" case, under the additional non- 
degeneracy assumption (1.9), it is known that any pulsating front with speed 
c is increasing in time and c > c*(e) (see [1]). 

The value c^A(j/(^) which appears in Theorem 1.6 is called the minimal 
speed of propagation of the pulsating travelling fronts propagating in the 
direction — e (satisfying the reaction-advection-diffusion problem (1.14)). 

We mention that the uniqueness of the pulsating travelling fronts, up to 
shifts in time, for each c > CQ^gj(e), has been proved recently by Hamel 
and Roques [7] for "KPP" nonlinearities. On the other hand, a variational 
formula for the minimal speed of propagation CQ^^j(e), in the case of a 
KPP nonlinearity, was proved in Berestycki, Hamel, Nadirashvili [2]. This 
formula involves eigenvalue problems and gives the value of the minimal 
speed in terms of the domain ^2 and in terms of the coefficients A, g, and / 
appearing in problem (1.14). The asymptotic behaviors and the variations 
of the minimal speed of propagation, as a function of the diffusion, advection 
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and reaction factors and as a function of the periodicity parameters, were 
widely studied in Berestycki, Hamel, Nadirashvili [3], El Smaily [4], El Smaily, 
Hamel, Roques [5], Heinze [8], Ryzhik, Zlatos [13], and Zlatos [20]. 



1.2 Main results 

In the periodic framework, having (in Theorems 1.5 and 1.6) the existence 
results and some qualitative properties of the pulsating travelling fronts prop- 
agating in the direction of a fixed unit vector — e G M'^, we search a variational 
formula for the unique speed of propagation c = c(e) whenever / is of the 
"combustion" type, and for the minimal speed c* = c^^^ j(e) whenever / is 
of the "ZFK" or the "KPP" type. We will answer the above investigations 
in the following theorem, but before this, we introduce the following 

Notation 1.7 For each function = (f){s,x,y) in C^'^(M. x Q) {for some 
6 G [0,1)), let 

F[(j)] := V^,y{AV^,y(j)) + {eAe)(f)ss + V^,yiAi(f)s)+ds{eAV^,y<j)) inV'{RxQ), 

where e = (e, 0, ■ ■ ■ , 0) G and e denotes a unit vector of M'^. 

The first main result deals with the "combustion" case. 

Theorem 1.8 Let e a unit vector o/M'^. Assume that Q is a domain satis- 
fying (1-2) and f is a nonlinear source satisfying (1-5) and (1-6). Assume 
furthermore that A and q satisfy (1-3) and (1-4) respectively. Consider the 
set of functions 

E = [if = Lp(s, x,y), if is of class C^'^(M x ^l) for each /i G [0, 1), 
F[(p] G C(M X f2), (y9 zs L— periodic with respect to x, (fg > 
m M X r2, Lp{—oo, ., .) = 0, (/3(+oo, ., .) = 1 uniformly in Q, and 
V ■ yl(Vx,j/V^ + ^v^s) =0 on M X dVL\ . 

For each Lp & we define the function R(p E C{Rx Q) as, for all (s, x, y) G 

R X n, 

D ( N F[Lp\{s,x,y)+ q-Va^^y^{s,x,y) + f{x,y,^) 
Rip{s,x,y)= — — ^ ^ + q{x,y)-e. 

If f is a nonlinearity of "combustion" type satisfying (1-7), then the unique 
speed c(e) that corresponds to problem (1-14) is given by 

c(e) = min sup Rip{s,x,y), (1-15) 
'^^^ (s,x,y)mxn 
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c(e) = max inf _Rip{s,x,y). (1-16) 

Furthermore, the min in (1-15) and the max in (1-16) are attained by, and 

only by, the function (j){s,x,y) = u (I'lT'd its shifts (j){s + T,x,y) 

for any r G M, where u is the solution of (l-M) with a speed c(e) (whose 
existence and uniqueness up to a translation in t follow from Theorem 1.5). 

The second result is concerned with "ZFK" nonhnearities. 

Theorem 1.9 Under the same notations of Theorem 1.8, if f is a nonlin- 
earity of "ZFK" type satisfying (1.5-1.6) and (1.8), then the minimal speed 
ch,A,q,f{e) given by 

= "^ig ^^P _ ^ 3;, y). (1.17) 

"^^ {s,x,y)mxn 



Furthermore, the min is attained by the function 0*(s, x, y) = u* (^^i^; 2:, y 

and its shifts (f)*{s + r, x, y) for any r G M, where u* is any solution of (1.1 4) 
propagating with the speed c*(e) = c^,^,qj(e). 

In particular, Theorem 1.9 yields that formula (1.17) holds in the "KPP" 
case (1.10) as well. 

Remark 1.10 In Theorem 1.8, the min and the max are attained by, and 
only by, the pulsating front (j){s, x, y) and its shifts (f){s + T, x, y) for all r G M. 
In Theorem 1.9, the min is achieved by the front (j)*{s,x,y) with the speed 
c*(e) and all its shifts 0*(s + T,x,y). Actually, if the pulsating front 0* is 
unique up to shift, then 0* and its shifts are the unique minimizers in formula 
(1.17). The uniqueness is known in the "KPP" case (see Hamel, Roques [7]), 
but it is still open in the general "ZFK" case. 

We mention that a max- min formula of the type (1.16) can not hold 
for the minimal speed c*(e) in the "ZFK" or the "KPP" case. A simple 
justification is given in section 2. 

The variational formulations of the speeds of propagation which are given 
in Theorems 1.8 and 1.9 are more general than those in Hamel [6] and Heinze, 
Papanicolaou, Stevens [9]. In Theorems 1.8 and 1.9, we consider nonhomo- 
geneous nonhnearities f = f{x, y, u) and the domain Vt is in the most general 
periodic situation. However, in [6] , the domain was an infinite cylinder of 
and the advection q was in the form of shear flows. Moreover, in this paper, 
the nonhomogeneous operator V ■ [AWu] replaces the Laplace operator Au 
taken in [(:>]. On the other hand, in [9], the domain VL was an infinite cylinder 
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in with a bounded cross section. Namely, ^2 = M x C where the 
cross section is a bounded domain in M^~^. Moreover, the authors did not 
consider an advection field in [9]. Finally, concerning the nonlinearities, they 
were depending only on u (i.e / = f{u) and is satisfying either (1.11) or 
(1.12)) in both of [6] and [9]. 

Besides the fact that we consider here a wider family of diffusion and 
reaction coefficients, our assumptions are less strict than those supposed in 
[9] and [16]. Roughly speaking, the authors, in [9] and [16], assume a stability 
condition on the pulsating travelling fronts. We mention that such a stability 
condition is fulfilled in the homogenous setting; however, it has not been 
rigorously proved so far that this condition is satisfied in the heterogenous 
setting. Meanwhile, the assumptions of the present paper only involve the 
coefficients of the reaction-advection-diffusion equation (1.14), and they can 
then be checked easily. 

Actually, in the "KPP" case, another "simpler" variational formula for 
the minimal speed c*(e) = ^ ^ j(e) is known. This known formula involves 
only the linearized nonlinearity / at m = 0. Namely, it follows from [2] that 

Theorem 1.11 (Berestycki, Hamel, Nadirashvili [2]) Let e be a fixed 
unit vector in and let e = (e, 0, . . . , 0) G M^. Assume that f is a "KPP" 
nonlinearity and that Q,A and q satisfy (1-2), (1-3) and (1-4) respectively. 
Then, the minimal speed c*(e) of pulsating fronts solving (1-14) o,nd propa- 
gating in the direction of e is given by 

c*(e) = c^.A,,,/(e) = min^, (1.18) 

' A>0 A 

where k{\) = kn^e,A,q,(W is the principal eigenvalue of the operator LQ^e,A,q,c,x 
which is defined by 

+ [\^eAe + \V-{Ae) + \q-e+(\^l) ^ ' 

acting on the set 

E\= {'?/'€ C'^{Vl),iI) is L-periodic with respect to x and 
V ■ A^SJil) = ■ Ae)i) on dVl} . 

In our last result, we prove that formula (1.17) implies formula (1.18) in 
the "KPP" case, but under some additional assumptions on the advection 
and the diffusion coefficients. This result gives an alternate proof of the 
formula (1.18). 
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Theorem 1.12 Let e be a fixed unit vector in M.'^ and let e = (e, 0, . . . , 0) G 
M^. Assume that f is a "KPP" nonlinearity and that A and q satisfy (1-2), 

(1.3) and (1-4) respectively. Assume, furthermore, that v ■ Ae = on dQ 
(in the case where dQ 7^ 0^'. Then, formula (1.17) implies formula (1.18). 

2 Main tools: change of variables and maxi- 
mum principles 

In this section, we introduce some tools that will be used in different places 
of this paper in order to prove the main results. 

Throughout this paper, e will denote the vector in defined by 

e = (e,0,---,0) = (e\---,e^O,---,0), 

where e^, ■ ■ ■ , e'^ are the components of the vector e. 

Our study is concerned with the model (1.14). Having a "combustion", 
a "ZFK", or a "KPP" nonlinearity, together with the assumptions (1.3) and 

(1.4) , problem (1.14) has at least a classical solution (c, u) such that c > and 
Ut > (see Theorems 1.5 and 1.6). The function u is globally C^''^(]R x Q) 
and C^'f^ with respect to (x, y) variables (for every fi G [0,1)). It follows 
that \/x,y{A'Vu) G C(M X r2). Having a unit direction e G M."^, and having a 
bounded classical solution (c, u) of (1.14) with c = c(e) (combustion case) or 
c > c*(e) (ZFK or KPP case), we make the same change of variables as Xin 
[19]. Namely, let (p = (l){s,x,y) be the function defined by 

(g — X • c \ 
,x,yj for all s G M and {x,y) G Q. (2.1) 

Then, for all (s, x, y) G M x f2, 

[Vx,j, ■ (AV^,y0) + ieAe)(l)ss + ^x,y ■ (^e0,) + dsieAW^^ycj)) ] (s, x, y) 

= V,,3, ■ {AVu){t,x,y), 

where s = x ■ e + ct. Consequently, 

F[0](s, X, y) = V^,y ■ (AV^,j,(/>) + {eAe)(f)ss + V^,^, ■ iAecf)^) + a,(eAV^,j,0) 

is defined at each point (s, x,y) eM.xQ and the map (s, x, y) t— > F[0](s, x, y) 
belongs to C(]R x Tt). 
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In all this paper, L = will denote the operator acting on the set E 
(given in Theorem 1.8) and which is defined by 

Llp = Vx-y ■ {AVx,yy:>) + {eAe)(fss + Vx,y ■ (Aeifs) + ds{eAVx,yy:>) 

+ q ■ ^x,yV + (<? ■ e - c)lPs in C(M x H) (2.2) 
= + q ■ Vx,y<^ + (q' ■ e — c)ips in C(M x Q), 

for all (f & E. 

It follows from above that if = (f){s,x,y) is a function that is given 
by a pulsating travelling {c,u) solving (1.14) (under the change of variables 
(2.1)), then F[<p] e C(M x H), is globally bounded in C^'^'{R x H) (for every 
fi G [0, 1)) and it satisfies the following degenerate elliptic equation 

L(j){s,x,y) + f{x,y,(f)) = F[(j)]{s,x,y) + q ■ Vx,y(j){s,x,y) 

+ {q- e - c)(f)sis,x,y) + f{x,y,(l)) =0 

in M X fi, together with the boundary and periodicity conditions 

J is L— periodic with respect to x, (o a\ 

\ V ■ A(Va:,j;0 + e0,) = on M X n. ^ ' ^ 

Moreover, since uit, x, y) as x-e ^ — cxd and n(t, a;, y) ^ 1 as x-e ^ +oo 
locally in t and uniformly in y and in the directions of M*^ which are orthogonal 
to e, and since is L— periodic with respect to x, the change of variables 
s = X ■ e + ct guarantees that 



0(— oo, ., .) = and 0(+oo, ., .) = 1 uniformly in (x, y) G fi. (2.5) 
Therefore, one can conclude that G -E. 

Remark 2.1 It is now clear that a max-min formula of the type (1-16) can 
not hold for the minimal speed c*(e) > m the "ZFK" or the "KPP" case. 
Indeed, for each speed c > c*(e), there is a solution {c,u) of (1.14) such that 
Ut > 0, which gives birth to a function = (f){s,x,y) under the change of 
variables (2.1). Owing to the above discussions the function cf) E E and it 
satisfies 

c = R(l){s, X, y) for all [s, x, y) G M x 

Therefore 

sup inf _Rip{s,x,y) > c. 

ipGE {s,x,y)eRxQ 
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Since one can choose any c > c*(e), one concludes that 

sup inf _R(p{s,x,y) = +00 

ip&E (s,a;,j/)eKxSl 

m the "ZFK" or the "KPP" case. 

Remark 2.2 (The same formulae, but over a subset of E) If the restric- 
tion of the nonlinear source f in (1-14) is C^'^iVL x [0, 1]), one can then con- 
clude that (see the proof of Proposition 6.3 in [1]) any solution u of (I.I4) 
satisfies: 

y{t,x,y) G M X n, \duu{t,x,y)\ < Mdtu{t,x,y) 
for some constant M independent of (t,x,y). In other words, the function 

X, y) = u{{s - X ■ e)/c, x, y) 

(where c = c(e) in the "combustion" case, and c = c*(e) in the "ZFK" or the 
"KPP" case) satisfies 

y{s,x,y) e M X n, \dss^{s,x,y)\ < (M/c) ds^{s,x,y). 

Let E' be the functional subset of E defined by 

E' = {^e E, 3C > 0, V(s,x,y) G M x H, x, y)| < C ds'^{s,x,y)} . 

The previous facts together with the discussions at the beginning of this sec- 
tion imply that the functions (j) o,nd (f)* of Theorems 1.8 and 1.9 are elements 
of E G E. These theorems also yield that the max-min and the min-max 
formulcB can also hold over the subset E of E. 

Namely, in the case of a "combustion" nonlinearity 

c(e) = min sup Rip{s,x,y) (2.6) 
vs-b' [s,x,y)mxn 

and 

c(e) = max inf _Rip{s,x,y). (2.7) 

ifiGE' {s,x,y)eRxn 

Moreover, the min and the max are attained at, and only at, the function 
(f){s, X, y) and its shifts (f){s + r, x, y) for any r G M. 

On the other hand, only a min-max formula holds in the case of "ZFK" 
or "KPP" nonlinearities. That is 

c*(e) = min sup Rip{s,x,y). (2.8) 

V^^' (s,x,y)GRxn 

Moreover, the min is attained at the function 4>*{s, x, y) and its shifts (f)*{s + 
r, X, y) for any r G M. 
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In the proofs of the variational formula which were given in Theorem 
1.8 and Theorem 1.9, we will use two versions of the maximum principle in 
unbounded domains for some problems related to (2.2-2.4) and (2.5). Such 
generalized maximum principles were proved in Berestycki, Hamel [1] in a 
slightly more general framework: 

Lemma 2.3 ([ ]) Let e be a fixed unit vector in Mf^. Let g{x, y, u) he a glob- 
ally hounded and globally Lip schitz- continuous function defined m ^2 x M and 
assume that g is non-increasing with respect to u in Q x (— oo, S] for some S > 
0. Let G M and := (— oo, h) x Q. Let c 7^ and (f)^{s, x, y), (f)'^{s, x, y) 

be two bounded and globally C^'^ ( j functions (for some fi > 0) such that 



L0i + (7(x,y,0i) > m P'(S-), 
L<P' + g{x,y,(P^) < zn P'(S^), 



■ A [e{(Pl - (Pi) + V^,j,((/>^ - 02)] < on (-00, h] x dn, (2-9) 
lim sup _ [0^(s,x,?/) - 02(s,a;,?/)] < 0, 

so~^ °° {s<sa, {x,y)en} 



where 



L(f) := V^,y{AV^,y(f)) + {eAi)^ss + '^x,y{Ai^s) + ds{eAV^,y(f)) 
+ q- Vx,y(j) + (g ■ e - c)(f)s 



■'S-i 

1>N 



and e denotes the vector (e, 0, ■ ■ ■ ,0) G 

If (f^ ^ ^ iiT' (ind (j)^{h,x,y) < (j)'^{h,x,y) for all {x,y) G then 

4>^ < 0^ in S^. 

Remark 2.4 Note here that (p^, 0^, q, A and g are not assumed to be L— periodic 
in X and that q is not assumed to satisfy (1.4)- 

Changing (j)^{s,x,y), (j)'^{s,x,y) and g{x,y,s) into 1 — (f)^{—s,x,y), 1 — 
(f)"^ {— s , X , y) and —g{x,y,l — s) respectively in Lemma 2.3 leads to the fol- 
lowing 

Lemma 2.5 ([ 1 ]) Let e be a fixed unit vector in Mf^. Let g{x, y, u) be a glob- 
ally bounded and globally Lip schitz- continuous function defined m x M and 
assume that g is non-increasing with respect to m m x [1 — 5, +00) for some 
(5 > 0. Let G M and := (/i, +00) x Vt. Let c 7^ and (j)^{s^ x, y), ^^(s, x, y) 
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be two bounded and globally C^'^ ( j functions (for some yU > such that 



L<f)' + g{x,y,<P') > zn 

i^-A[e{(j)l-(f)^,)+V,,y{(l)^-(p^)] < on [/i,+cx)) X 9fi(2- 11) 
lim sup _ - < 0, 

.0^+00 ^^^^^^^ {x,y)(^n} 

where L is the same operator as in Lemma 2.3. 

If (f)"^ > 1 ~ S in and (j)^{h, x, y) < x, y) for all {x, y) G Q, then 

0^ < 0^ m S+. 



3 Case of a "combustion" nonlinearity 

This section is devoted to prove Theorem 1.8, where the nonhnearity / sat- 
isfies the assumptions (1.5-1.6) and (1.7). 

3.1 Proof of formula (1.15) 

Having a prefixed unit direction e G M"^, and since the coefficients A and q 
of problem (1.14) satisfy the assumptions (1.3) and (1.4), it follows, from 
Theorem 1.5, that there exists a unique pulsating travelling front (c(e),u) 
{u is unique up to a translation in the time variable) which solves problem 
(1.14). Moreover, dfU > in M x f2. We will complete the proof of (1.15) via 
two steps. 

Step 1. After the discussions done in the section 2, the existence of a 
classical solution (c(e), m), satisfying (1.14), implies the existence of a globally 
C^{M. X Vt) function (j){s,x,y) satisfying O<0<linMxn, with 

4> is L— periodic with respect to x, 
L<j){s, X, y) + /(x, y, 0) = in P'(M x H), 

V ■ A{Vx^y(p + ecps) = in M X ^ ' 

0(— oo, ., .) = 0, and 0(+oo, ., .) = 1 uniformly in (x, y) G fi, 

where L is the operator defined in (2.2) for c = c(e). We also recall that the 
two functions u and satisfy the relation 

u{t, X, y) = (f>{x ■ e + c{e)t, x, y), (t, x,y) & M. x Q. 

One has ds4> > in M x 1] and this is equivalent to say that the function 
u = u{t,x,y) is increasing in t, since c(e) > 0. 
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Together with the facts in section 2.1, one gets that the function (f) G E. 
Furthermore, (3.1) yields that 

Vs G M, V(x,i/) e n, c(e) = R(f){s,x,y), (3.2) 

and 

L0(s,x,?/) + /(x,?/,0) = 0, (3.3) 

where R(j) is the function defined in Theorem 1.8. In other words, the 
L— periodic (with respect to x) function R(j) is constant over M x f2 and 
it is equal to c(e). 

It follows, from (3.2) and from the above explanations, that 

c(e) > inf sup Rip{s,x,y). 

"P^^ (s,a;,j/)eMxf! 

To complete the proof of formula (1.15), we assume that 

c(e) > inf sup Rip{s,x,y). 

'^^^ (s,x,y)mxn 

Then, there exists a function ip = ip{s, x,y) E E such that 

c(e) > sup Rip{s,x,y). 

{s,x,y)eRxQ 

Since the function ip E E, one then has ^l>s{sy 2;, y) > for all (s, x,y) eM.xQ. 
This yields that 

Lilj{s,x,y) + fix, y,ilj) <0 inR xTi, (3.4) 

where L is the operator defined in (2.2) for c = c(e). 

Notice that the later holds for each function of the type 

^^{s,x,y) := tpis + T,x,y) 

because of the invariance of (3.4) with respect to s and because the advection 
field q and the diffusion matrix A depend on the variables {x, y) only. That 
is 

L^p^is, X, y) + /(x, y, ?^^) < in M x H. (3.5) 

Step 2. In order to draw a contradiction, we are going to slide the function 
ip with respect to (p. From the limiting conditions satisfied by these two 
functions, there exists a real number i? > such that 

J (p{s,x,y) < 9 for all s < —B, {x,y) G fl, 
\ ip{s, X, y) > 1 — p for all s > B, (x, y) G f2, 
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and 



4>{B, x,y) > 1 — p for all (x, y) G Q, 



(3.6) 



where 6 and p are the values that appear in the conditions (1.7) satisfied by 
the "combustion" nonlinearity /. Taking r > 2B, and since is increasing 
with respect to s, one gets that (f){—B,x,y) < tp'^{—B,x,y) for all {x,y) G Q 
and ■ip'^ > 1 — p in Sl^. 

It follows from Lemma 2.3 (take 6 = 6, h = —B, cfr" = (p, and (j)"^ = tp'^) 
that (p < ip'^ in Sl^. Moreover, Lemma 2.5 (take 6 = p, h = —B, cp^ = 
0, and (p'^ = ip'^ ) implies that (p < ip'^ in Sl^. Consequently, (p < ip'^ in 
M X n for all r > 2B. 

Let us now decrease r and set 



First one notes that r* < 2B. On the other hand, the limiting conditions 
ip{—oo,.,.) = and 0(+oo,.,.) = 1 imply that r* is finite. By continu- 
ity, <P < in M X f2. Two cases may occur according to the value of 



Since the functions ip and (p are globally C^(]R x Q) there exists r] > such 
that the above inequality holds for all r G [r* — rj, t*]. Choosing any r in the 
interval [r* — rj, r*], and applying Lemma 2.3 to the functions ip'^ and (p, one 
gets that 

(p{s,x,y) < ip'^{s,x,y) for all s < —B, {x,y) G fl, 

together with the inequality 

(p{s, X, y) < ip'^{s, X, y) for all s G [—B, B], and for all (x, y) G fl. 

Owing to (3.6) and to the above inequality, it follows that 

ip^{B, x,y) > 1 — p in Q. 

Moreover, since the function ip is increasing in s, one gets that ip^ > 1 — p 
in S^. Lemma 2.5, applied to (p and ip'^ in Ej, yields that 



T* = inf{r G M, (p < in R x Q }. 



sup [(p — W ) ■ 

[-B,B]xn 

case 1: suppose that 




sup 

[~B,B]xTl 




(p{s, X, y) < ip^i^s, X, y) for all s > B, (x, y) G Q. 
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As a consequence, one obtains < in M x ^2, and that contradicts the 
minimahty of r*. Therefore, case 1 is ruled out. 
case 2: suppose that 

sup (0-^^*)=O. 

[~B,B]xTi 

Then, there exists a sequence of points {sn,Xn,yn) in [—B,B] x Q such that 

(j){Sn, Xn, Vn) " V {Sn, Xn, ?/„) ^ aS n +00. 

Due to the L— periodicity of the functions </> and tp, one can assume that 
{xn, Vn) G C. Consequently, one can assume, up to extraction of a subse- 
quence, that Xn, Vn) {s, X, y) G [— -B, -B] x C as n — >• +00. By continu- 
ity, one gets (pis, x, y) = ^Ij^*{s, x, y). 

We return now to the variables {t,x,y). Let 

z{t,x,y) = 4){x ■ e + c{e)t,x,y) — ip{x ■ e + c{e)t + T* ,x,y) 

= u{t, X, y) — ip{x ■ e + c(e) t + r*, x, y) for all (t, x, y) G M x fi. 

Since the functions cf) and ip are in it follows that the function z is globally 
C^{R X Vt) and it satisfies 

V(t,x,y) gMx fi, V.,^y{AVz){t,x,y) = F[(P]{s,x,y)-F[iP^'']{s,x,y), 

where s = x-e + c{e)t. Thus, Vx,y- {AVz) e C(Mxf2). Moreover, the function 
z is non positive and it vanishes at the point {{s — x ■ e) / c{e) , x, y) . It satisfies 
the boundary condition u ■ (AVz) = on M x dfl. Furthermore, it follows, 
from (3.2) and (3.4), that 

dtz - V^,y ■ (AVz) + q{x, y) ■ V^^yZ < f{x, y, 0) - f{x, y, ip'"'). 

However, the function / is globally Lipschitz-continuous in f2 x M; hence, 
there exists a bounded function h{t, x, y) such that 

dtZ — Vx,y ■ (AVz) + q{x, y) ■ Vx,yZ + b{t, x,y) z < in R x Q, 

with z{t, x,y) < for all (t, x, y) G M x fi. 

Applying the strong parabolic maximum principle and Hopf lemma, one 
gets that z{t,x,y) = for all t < {s — x-e)/c{e) and for all {x,y) G VL. On the 
other hand, it follows from the definition of z and from the L— periodicity of 
the functions and ip that ^(t, x,y) =0 for all (t, x,y) G Mx Consequently, 

0(s, X, y) = ip^ (s, X, y) = ip{s + r*, x, y) for all (s, x, y) G M x 1]. 

Referring to the equations (3.3) and (3.5), one gets a contradiction. Thus, 
case 2 is ruled out too, and that completes the proof of the formula (1.15). 
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Remark 3.1 (The uniqueness, up to a shift, of the minimizer in (1.15)) 

If il) E E is a minimizer in (1.15). The above arguments imply that case 2 
necessarily occurs, and that ip is equal to a shift of 0. In other words, the 
minimum in (1.15) is realized by and only by the shifts of (p. 

3.2 Proof of formula (1.16) 

In this subsection, we are going to prove the "max-min" formula of the speed 
of propagation c(e) whenever the nonlinearity / is of the "combustion" type. 
The tools and techniques which one uses here are similar to those used in the 
previous subsection. However, we are going to sketch the proof of formula 
(1.16) for the sake of completeness. 

As it was justified in the previous subsection, one easily gets that 

c(e) < sup inf _R(p{s,x,y) 

and 

V (s, X, y) G M X r2, c(e) = R(f>{s, x, y), 

where 

0(s, x,y) = u{{s — X ■ e)/c(e), X, y) , for all (s, x, y) G M x fi, 

and u = u{t, x, y) is the unique (up to a translation in t) pulsating travelling 
front solving problem (1.14) and propagating in the speed c(e). We recall that 
the function (j) E E (see section 2). It follows that the function satisfies the 
following 

is L— periodic with respect to x, 
is increasing in s G M, 
< L0(s,x,|/) + /(x,|/,0) = in R X n, (3.7) 
u ■ A(Va;,j/0 + e0j = in M X 

0(— oo, ., .) = 0, and 0(+oo, ., .) = 1 uniformly in {x,y) G Q, 

where L is the operator defined in (2.2) for c = c(e). 

Notice that the later holds also for each function of the type 

0^(s,a;,?/) := 0(s + T,x,y) 

because of the invariance of (3.8) with respect to s and because the advection 
field q and the diffusion matrix A depend on the variables (x, y) only. 
To complete the proof of formula (1.16), we assume that 

c(e) < sup inf _Rip{s,x,y). 
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Hence, there exists G E such that 

c(e) < Rip (s, X, y), for all (s, x,y) E M. x Q. 

Since the function tp E E, one then has ipsi^: 2;, y) > for all (s, x, y) G Mx 
This yields that 

Lip{s, X, y) + f{x,y,tp) > inRxQ. (3.8) 

To get a contradiction, we are going to slide the function (p with respect 
to tp. In fact, the limiting conditions satisfied by ip and (p, which are elements 
of E, yield that there exists a real positive number B such that 

J ip{s,x,y) < 6 for all s < —B, G fi, 

\ X, y) > 1 — p for all s > B, (x, y) G fi, 

and 

ij{B, x,y)>l-p for all (x, y) G H, (3.9) 

where 6 and p are the values appearing in the conditions (1.7) satisfied by the 
nonlinearity /. Having r > 2B, one applies Lemma 2.3 (taking 6 = 6, h = 
—B, (p^ = tp, and 0^ = 0^) and Lemma 2.5 (taking 6 = p, h = —B, (p^ = 
ip, and 0^ = 0"^) to the functions (p'^ and over the domains Sl^ and 
respectively, to get that tp < (p'^ in Sl^ and < (^"^ in S^^. Consequently, 
one can conclude that 

Vr >2B, tp <(p^ in M X H. 

Let us now decrease r and set 

T* = inf{r G M, < 0^ in M x H }. 

It follows, from the limiting conditions ip{+oo, ., .) = 1 and (p{—oo, ., .) = 0, 
that r* is finite. By continuity, we have ip < (p'^* . In this situation, two cases 
may occur. Namely, 

case A: sup — 0'^ )< 0, 

[~B,B]yTi 

or 

case B: sup (^/; — 0"^ ) = 0. 

[-B,B]yH 

Imitating the ideas and the skills used in case 1 and case 2 during the proof 
of formula (1.15), one gets that case A (owing to minimality of r*) and case 
B (owing to (3.7) and (3.8)) are ruled out. 
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Therefore, the assumption that 

c(e) < sup inf _Rip{s,x,y) 

is false, and that completes the proof of formula (1.16). 

Remark 3.2 (The uniqueness, up to a shift, of the maximizer in (1. 

Similar to what we have already mentioned in Remark 3.1, if tp & E is a 
maximizer in (1.16), then the above arguments yield that case B necessar- 
ily occurs, and that ip is equal to a shift of (p. One then concludes that the 
maximum in (1.16) is realized by, and only by, the shifts of (p. 



4 Case of "ZFK" or "KPP" nonlinearities: 
proof of formula (1.17) 

This section is devoted to the proof of Theorem 1.9. We assume that the 
nonlinear source / is of "ZFK" type. Remember that this case includes 
the class of "KPP" nonlinearities. Namely, / = f{x, y, u) is a nonlinearity 
satisfying (1.5-1.6) and (1.8). We will divide the proof of formula (1.17) into 
3 steps: 

Step 1. Under the assumptions (1.2), (1.3), and (1.4) on the domain f2, 
the diffusion matrix A, and the advection field q respectively, and having 
a nonlinearity / satisfying the above assumptions. Theorem 1.6 yields that 
for c = CQ^,jj(e), there exists a solution u* = u*(t,x,y) of (1.14) such that 
u^it, x,y) > for all (t, x,y) G M x fi. In other words, the function 0* defined 

by 

— X, yj , {s,x,y) eRxQ 
is increasing in s G M. Owing to section 2, 0* satisfies 

F[4>*] + q ■ V.>* + {q-e - c*(e))0:, +/(x, y,4>*) = in M x H (4.1) 
together with boundary and periodicity conditions 



0* is L— periodic with respect to x, 
V ■ A(Vx.,j;0* + e0*) = on R X n. 



(4.2) 



Moreover, (4.1) implies that 

V (s,x,?/) G M X n. 



^ F 0* s, X, y + g ■ V^,^0* g, X, y) + f{x, y, 0* 

^i^)= jnj? ^ + Q{x,y) 

ds(j)*{s,x,y) 

= R(f)*{s,x,y), 
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and hence 



c*(e)>inf sup F[^]i^^^^y)+<l-^^,y^ + fi^^y^v) ^q^^^yy,_ 
^e^(.,x,s;)GRxn ds(f){s,x,y) 

In order to prove equality, we argue by contradiction. Assuming that the 
above inequahty is strict, one can find 6 > such that 

c*(e)-5>inf sup mis,x,y)+ g.V.,y^ + f{x,y,^) 
^eS(.,x,y)eRxT7 dsip{s,x,y) 

(4.4) 

To draw a contradiction, we are going to approach the "ZFK" nonhnearity / 
by a sequence of "combustion" nonhnearities (/e)^ and the minimal speed of 
propagation by the sequence of speeds {cg)o corresponding to the functions 
{fg)e- The details will appear in the next step. 

Step 2. Let x be a C^(M) function such that < x < 1 in M, x(^) = 
for all M < 1, < x{u) < 1 for all u G (1,2) and x(«) = 1 for all u > 2. 
Assume moreover that x is non-decreasing in R. For all 9 G (0, 1/2), let xe 
be the function defined by 

VmGM, xe{u) = x{u/e). 

The function xe is such that < < 1, < xe < 1 in {—oo, 6^], < xe < 1 
in (6,26) and xe = 1 in [2^, +oo). Furthermore, the functions xe are non- 
increasing with respect to 9, namely, xoi ^ Xe2 if < 6'i < 6^2 < 1/2. 
We set 

fe{x, y, u) = f{x, y, u) xeiu) for all (x, y,u) eilxR. 

In other words, we cut off the source term / near u = 0. 

For each 9 G (0, 1/2), the function /g is a nonlinearity of "combustion" 
tj^e that satisfies (1.5-1.6) and (1.7) with the ignition temperature 9. There- 
fore, Theorem 1.5 yields that the existence of a classical solution {cg,ug) of 
(1.14) with the nonlinearity fg. Furthermore, the function ue is increasing in 
t and unique up to translation in t and the speed cg is unique and positive. 

It was proved, through Lemma 6.1 and Lemma 6.2 in Berestycki, Hamel 
[1], that the speeds cg are non-increasing with respect to 9 and 

Cg / c*(e) as 6^ \ 0. 

Consider a sequence 6'„ \ 0. Then, there exists rio G N such that cg^ > 
c*(e) — S for all n > uq (or equivalently 9n < 6'„p). 
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In what follows, we fix 6 such that 6 < 6^^- One consequently gets cg > 
c*(e) — 6. On the other hand, it follows, from the construction of fg, that 
/ > /e in X R. Together with (4.4), one obtains 

cg > mf sup — ^-^ — h g ■ e. (4.5) 

Thus, there exists a function ip & E such that 

F[ip]{s,x,y) + q-V^^y-^is.x.y) + fg{x,y,ilj) 
cg> — 7^ ^ + q{x,y)-e. 4.6 

However, ifjg^s.x.y) > for all {s,x,y) G M x f2. Thus, the inequality (4.6) 
can be rewritten as 

L^(s, X, y) + fe{x, y,^p)<0 inRxn, (4.7) 

with ip E E and L is the operator defined in (2.2) for c = cg. 
For each r G M, we define the function ip'^ by 

ip^{s, X, y) = ip{s + T, X, y) for all (s, x, y) G M x 

Since the coefficients of L are independent of s, the later inequality also holds 
for all functions ip'^ with r G M. That is, 

W{s, X, y) + fe{x, y, r)<^ m M x H. (4.8) 



Step 3. For the fixed 9 (in step 2), the function fg is a "combustion" 
nonlinearity whose ignition temperature is 9. There corresponds a solution 
{cg,ug) of (1.14) within the nonlinear source fg. We define (pg by 

, x,y], for all (s, x, y) G M x 

Cg J 

Referring to section 2, one knows that (pg E E and thus it satisfies the 
following equation 

L(Pg{s,x,y) + fe{x,y,(Pg) = in M x H. (4.9) 

Now, the situation is exactly the same as that in step 2 of the proof of 
formula (1.15) because the nonlinearity fg is of "combustion" type. The little 
difference is that / (in step 2 of the proof of formula (1.15)) is replaced here 
by fg, and the function of equation (3.3) is replaced by the function (pg of 
(4.9). Thus, following the arguments of subsection 3.1 and using the same 
tools of "step 2" as in the proof of formula (1.15), one gets that the (4.4) is 
impossible and that completes the proof of formula (1.17). □ 
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Remark 4.1 We found that one can use another argument (details are be- 
low) different from the sliding method in order to prove the min — max for- 
mul(B for the speeds of propagation whenever f is a homogenous (i.e f = 
f{u)) nonlinearity of "combustion" or "ZFK"type and Q = M^. Meanwhile, 
the sliding method, that we used in the proofs of formulcB (1-15) and (1.17), 
is a unified argument that works in the general heterogenous periodic frame- 
work. 

Another proof of formulae (1.15) and (1.17) in a particular frame- 
work: 

Here, we assume that / = /(u), and Vt = M^. Following the same procedure 
of "step 1" in the previous proof, one gets the inequality 

c*(e) > inf sup Rip{s,x,y). 

Now, to prove the other sense of inequality, we assume that 
c*(e) > inf sup R(f{s,x,y), 

and we assume that / is of "ZFK" type ^. Then, as it was explained in "step 
2" of the previous proof, one can find ip ^ E, 6 > 0, 9 > 0, and d > such 
that 

c*(e) - 5 <d <C0 < c*{e) 

where 

V(s, x,y) eR xQ, d> c*{e) - 6 > R^p{s, x, y), 

and fe{u) = f{u) Xeiu) < f{u) for all m G M is of "combustion" type (cg is 
the speed of propagation, in the direction of — e, of pulsating travelling fronts 
solving (1.14) with the nonlinearity fe and the domain VL = M^). 
Hence, for all {t,x,y) G M x M^, 

d> ""[^J^^' ^' ; 7-^^^.^' + ^^^^^ + y) ■ e. (4.10) 

Let u{t, X, y) = ip{x-e-\-dt, x, y). As it was explained in section 2, the function 
u satisfies 

^ Ut-W ■ {A{x, y)Wu) - q{x, y)-Vu- fe{u) > 0, t G M, (x, y) G H, 
V ■ A Vu{t, x,y) = 0, t e M, (x, y) G dQ, 

A - k-e (4.11) 

VA; G ||LiZ, y{t,x,y) G M x 1], u{t ——,x,y) = u{t,x + k,y), 

(Jj 

1=1 

< u< 1. 



^The case where / is of "combustion" type follows in a similar way. 
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Let < UQ{x,y) < 1 be a function in C(M^) such that Uo{x,y) ^ as 
x-e —CO, and Uo{x, y) ^ 1 as x-e ^ +oo, uniformly in y and all directions 
of M'^ which are orthogonal to e. Let m be a pulsating front propagating in 
the direction of — e with the speed cg and solving the initial data problem 

Ut = V- iA{x, y)Vu) + y)-Vu+ fe{u), t > 0, (x, y) E H, 
uiO,x,y) = Uo{x,y), (4.12) 
u ■ A S/u{t, X, = 0, t e M, (x, y) e 90. 

Having fe{u) as a "combustion" nonlinearity, it follows from J. Xin [18] 
(Theorem 3.5) and Weinberger [17], that 

Vr > 0, lim sup n(t, x — cte, y) = uniformly in y, for every c > cg, 

*^+°° \x\<r 

and lim inf u{t, x — cte, y) = 1 uniformly in y, for every c < cq. 

t^+oo \x\<r 

(4.13) 

This means that the speed of propagation cq corresponding to (1-14) is 
equal to the spreading speed in the direction of — e when the nonlinearity is 
of "combustion" type and the initial data uq satisfies the above conditions. 

For all (t,x,y) G [0, +oo) x Q, let w(t,x,y) = u{t,x,y) — u{t,x,y). It 
follows, from (4.11) and (4.12), that 

Wt-W ■ {A{x, y)'Vw) - q{x, y) ■ Vw + bw > 0, t > 0, (x, y) G O, 
V(x,y) gH, wiO,x,y)>0, (4.14) 
u ■ A Ww{t, x,y) = 0, t e M, (x, y) G dil, 

for some b = b{t, x,y) G C(M x O). The parabolic maximum principle implies 
that w > in [0, +oo) x Q. In other words, 

V(t, X, y) G [0, +oo) X Q, u(t, x, y) < u(t, x, y). 

However, for all c> d, 

lim u{t, X — cte, y) = lim ip{x ■ e + {d — c)t, x — cte, y) = 

t— >+oo t— ►+00 

locally in x and uniformly in y (since ip E E). Consequently, 

Vc > d, Vr > 0, lim sup u(t, x — cte, y) = uniformly in y. 

Referring to (4.13), one concludes that d > cg which is impossible {d < cg). 
Therefore, our assumption that c*(e) > inf sup Rip{s,x,y) is false and 

that completes the proof of formula (1.17) in the case where / = f{u) and 
O = R^. □ 
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